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ABSTRACT 

We consider the time-dependent evolution of a relativistic jet following its breakout through the surface of a 
massive compact star, as envisaged in the collapsar model of gamma-ray bursts. At breakout, the jet is tightly 
collimated by the pressure of its hot cocoon, which is created as the jet traverses the star. After breakout, the 
cocoon pressure drops and the jet expands toward its natural opening angle. We show that the evolving opening 
angle of the jet produces a stratification of the total energy with the off-axis angle, resulting in a universal 
morphology. The angular structure is largely independent of the initial beam pattern and depends only on the 
luminosity of the central engine. With the minimal assumption of a constant energy release we reproduce the 
8~ 2 profile required to explain observations of afterglows. 

Subject headings: gamma-rays: bursts — hydrodynamics — supernovae: general 



1. INTRODUCTION 

After the observation of GRB 990123 and its afterglow, it 
was realized that Gamma-Ray Bursts (GRBs) are collimated 
and not isotropic explosions. GRB 990123 was a most prob- 
lematic event since it called for more than a solar mass of 
energy entirely converted into 7-ray photons (Kulkarni et al. 
1999) to explain its fluence under the assumption of isotropic 
emission. It also had a peculiar afterglow, with a steepening of 
the decay rate approximately one day after explosion (Castro- 
Tirado et al. 1999). This behavior had been predicted to be 
the signature of a collimated outflow (Rhoads 1997, 1999). 
After correcting for beaming the implied energy release was 
reduced to - 5 x 10 51 erg (Bloom, Frail & Kulkarni 2003), 
well below the energy crisis level of 10 53 erg. 

GRB jets were initially postulated to be uniform, i.e., with 
constant properties (energy flux, Lorentz factor, pressure) 
within a well-defined opening angle 6j and nothing outside. 
Such jets would produce a steepening afterglow light curve, 

with the break time scaling as 9j {E\ so /n) 1 ^ (Sari, Piran & 
Halpern 1999), where E[ so is the isotropic equivalent energy 
of the outflow and n the density of the ambient medium. Frail 
et al. (2001) and Panaitescu & Kumar (2001) discovered a 
remarkable anti-correlation between the jet break time and 
the isotropic equivalent energy release E 1 ^ in the prompt 
phase. Such a correlation, initially entirely empirical, can be 
accounted for if it is assumed that the same total energy is 
given to every GRB but is channeled into jets with different 
opening angles. A narrow jet would appear as a very bright 
GRB followed by an early-breaking afterglow, while a wide 
jet would produce a weak GRB and an afterglow with a late 
break. 

Rossi, Lazzati & Rees (2002) realized that the relation 
E 1 ^ S 2 j =constant may reflect a universal angular distribu- 
tion of jet energy rather than a distribution of opening angles 
among different jets (see also Lipunov, Postnov & Prokhorov 
2001; Zhang & Meszaros 2002). In this case one would write 

Rossi et al. also showed that such an energy pattern would 



produce an afterglow that reproduces the Frail et al. corre- 
lation and does not violate other observations. Different ob- 
served properties would reflect different observing geometries 
rather than different intrinsic jet properties. Several other jet 
morphologies have since been proposed and studied, specifi- 
cally general power-law profiles and Gaussian profiles (Gra- 
not & Kumar 2003; Salmonson 2003, Rossi et al. 2004; 
Zhang et al. 2004a). Despite this theoretical interest, research 
on structured jets has focused on their late-time dynamics 
(Kumar & Granot 2003) and observable properties (Perna, 
Sari & Frail 2003; Nakar, Granot & Guetta 2004) and not on 
the origin of the energy distribution. The only exception is the 
electromagnetic force-free model, in which a structured jet is 
predicted naturally (Lyutikov, Pariev & Blandford 2003). 

Some degree of structure in GRB jets is implied by their 
connection to supernova (SN) explosions (Galama et al. 1998; 
Stanek et al. 2003; Hjorth et al. 2003; Malesani et al. 2004). 
The shear forces and mixing of the cold stellar material with 
the outflowing relativistic plasma are expected to create an 
interface or hollow slower jet around the hyper-relativistic 
flow (MacFadyen & Woosley 1999; MacFadyen et al. 2001; 
Proga et al. 2003). The resulting energy distributions of the 
highly relativistic material are, however, far from a power-law 
(Zhang, Woosley & Heger 2004b). If one includes the trans- 
relativistic outflow energy a power-law profile is obtained, but 
with the steep index -3. 

In this letter we propose a mechanism to create structured 
jets with dE/aTl oc 8~ 2 as a result of the time evolution of the 
opening angle of an hydrodynamic jet. This mechanism re- 
lies on the nature of the jet propagation through the star but is 
more robust and predictable than the poorly understood shear 
effects and mixing instabilities. We show that, under a simple 
ansatz (constant luminosity engine) this mechanism yields a 
beam pattern following the prescription of eq. [2 under wide 
conditions. This letter is organized as follows: in § 2 we de- 
scribe the star, cocoon and jet conditions at breakout; in § 3 we 
compute the jet evolution during the release of the cocoon ma- 
terial and the resulting jet structure at large radii. We discuss 
our results and their implications in § 4 and we summarize 
our findings in § 5. 
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2. THE JET AND ITS COCOON AT BREAKOUT 

Consider a high entropy mildly relativistic jet (To > 1) gen- 
erated at some small radius ro ~ 10 7 — 10 s cm in the core of 
a massive star with an opening angle 9o- The jet tries to ac- 
celerate under the pressure of internal forces. We assume that 
the initial jet has enough momentum to define a direction of 
propagation and avoid a spherical explosion of the star. Such 
initial conditions are analogous to those adopted in numerical 
simulations (Zhang et al. 2003, 2004b). In the absence of any 
external material, it would accelerate in a broad conical out- 
flow, with its Lorentz factor scaling with radius as T oc r/ro. 

If the jet is surrounded by dense cold matter, as in the 
collapsar/hypernova scenario (Woosley 1993; MacFadyen & 
Woosley 1999), such acceleration is inhibited. Once the jet 
reaches supersonic speeds with respect to the stellar material, 
the ram pressure of the shocked gas ahead of the jet drives 
a reverse shock into the head of the jet, slowing the jet head 
down to subrelativistic speeds. The mixture of shocked jet 
and shocked stellar material surrounds the advancing jet to 
form a high pressure cocoon, analogous to the cocoon formed 
by a low-density jet from a radio galaxy advancing into the 
intergalactic medium (Scheuer 1974). As a result the jet dy- 
namics is governed by three competing effects. First, the in- 
ternal pressure accelerates the jet material, whose Lorentz fac- 
tor scales as V oc (E ; /Eo) I//2 if the acceleration is isentropic, 
where Ey is the jet cross section of a single jet and Eq its value 
at ro (see, e.g., Beloborodov 2003, § 3.1). If the jet suffers 
internal dissipation as it propagates, e.g., from entrainment of 
cocoon material or internal shocks, then the increase of V with 
E ; - is slower. This can be caused, e.g., by shear instabilities at 
the jet-cocoon boundary. Second, the head of the jet is slowed 
down by its interaction with the massive star. The larger the jet 
cross-section, the slower its head propagates. Both shocked 
jet and stellar material flow to the sides feeding the cocoon. 
Third, the cocoon pressure acts as a collimating force on the 
jet, which therefore becomes narrower and more penetrating, 
albeit less relativistic. 

Matzner (2003) developed a simple analytic treatment to 
follow this propagation. Using the Kompaneets approxima- 
tion (Kompaneets 1960) to describe the cocoon expansion, he 
was able to compute the jet propagation time within the star 
and the cocoon properties at the moment the jet head reaches 
the stellar surface — the breakout time, t\, T . For our purposes 
we need merely to compare the width of the cocoon at break- 
out with the width of the jet. 

Consider a jet with luminosity Lj that reaches the stellar 
surface at fb r . The energy stored in the cocoon (neglecting 
adiabatic losses) is E c = Ljit^-r^/c), where is the radius 
of the star. About half of the cocoon energy is transferred 
to the nonrelativistic shocked stellar plasma via the cocoon 
shock, with the other half remaining in relativistic, shocked 
jet material. If we adopt a relativistic equation of state for the 
entire cocoon, we can write the cocoon pressure as 

Pc ~ — (2) 
1 3V C ' 

where V c w r+ r\ is the volume and r j_ = v s h fbr the transverse 
radius of the cocoon (we assume that the jet occupies a neg- 
ligible fraction of the cocoon volume.) The velocity of the 
shock v s h driven by the cocoon pressure into the star can be 
computed by balancing the cocoon pressure against the ram 
pressure exerted by the stellar material (which has negligible 
internal pressure): v s h = J p c /p*, where p+ is the matter den- 



to obtain an equation for the cocoon pressure that reads (see 
also Begelman and Cioffi 1986 for a similar treatment in the 
context of extragalactic radio sources): 

/ . x 1 V 2 / j . „ \ 1/2 



Pc = 



L JP* 



3r*r br 



(3) 



We assume that the head of the jet propagates subrelativis- 
tically through the star and write fbr = T] r */ C < where rj > 1. 
Then, adopting the notation Q = WQ. K and using cgs units 
throughout, we have 
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The opening angle of the cocoon is given by 
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i.e., of the order of tens of degrees. 

Now let us consider the properties of the jet as it reaches 
the breakout radius. The jet pressure is given by pj = 
Lj/{AcYijY 2 j), where Tj is the Lorentz factor of the jet material 
before it is shocked at the jet head. Setting pj = p c (Begelman 
& Cioffi 1989; Kaiser & Alexnder 1997) with E; = nr^Oj, we 
find 
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For isentropic flow between ro and r* the confined jet reaches 
a Lorentz factor 



with somewhat lower values if the jet is dissipative. The cor- 
responding jet opening angle at breakout is < 1°; however, 
the fact that #/,b r Tj.br < 1 suggests that the jet should freely 
expand to a few times 0y.br after exiting the star. Even so, the 
jet at breakout is expected to be more than an order of mag- 
nitude narrower than the cocoon. A more detailed treatment 
of the evolution of the cocoon and its interaction with the jet 
can be obtained by numerical integration of the set of differ- 
ential equations that govern the flow. This allows us to take 
into account the density gradient of the star and the shape of 
the cocoon. Results (Lazzati et al. 2005) do not differ sub- 
stantially from these analytic estimates. 

These results may be altered if substantial dissipation due 
to recollimation shocks and/or shear instabilities takes place. 
We do not attempt to model these effects here. Recollimation 
shocks are however relatively weak, being mostly oblique to 
the flow. As a matter of fact, these results are in qualitative 
agreement with results from numerical simulations (Zhang, 
Woosley & MacFadyen 2003). Simulations show, indeed, that 
isentropic conditions are not respected throughout the propa- 
gation, since collimation shocks take place. These shocks ef- 
fectively create a new nozzle at a larger radius r > ro which 
modifies the scaling of the Lorentz factor with the jet cross 
section. What is of relevance here is that no matter how wide 
the jet is initially (and it is likely to be poorly collimated, espe- 
cially if initial collimation is provided by the accretion disk), 
the jet that emerges from the star is narrow and highly col- 
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plying an exponentially increasing opening angle of the form 
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FIG. 1 . — Energy distribution for the afterglow phase for three instantaneous 
beam patterns (see inset). In all three cases a well defined dE jdQ. oc 0~ 2 sec- 
tion is clearly visible. Only the edge of the jet and its core show marginal 
differences. The results shown are for a jet/star with L = 10 s ' erg s~', 
?grb = 40 s and r* = 10" cm. Inset: Instantaneous beam patterns that reach 
the surface of the star. The solid line shows a uniform jet, dashed line shows 
a Gaussian energy distribution, while the dash-dotted line shows an edge 
brightened (or hollow) jet. 



all numerical simulations of jets in collapsars (MacFadyen & 
Woosley 1999; MacFadyen, Woosley & Heger 2001; Zhang 
et al. 2003). The cocoon, on the other hand, covers a wide 
solid angle, several tens of degrees across. In the next section, 
we explore what happens after jet breakout. 

3. JET BREAKOUT 

We now consider the jet development after the breakout, 
with the inner engine still active. At this stage, which has not 
been investigated in numerical simulations so far, it is likely 
that dissipation will have a lesser role, since the jet, as we 
shall see, is de-collimated rather than recollimated. 

As the jet reaches the stellar surface, it clears a channel 
for the cocoon. The cocoon material is therefore now free 
to expand out of the star and its pressure drops. We assume 
that from this moment on the shock between the cocoon and 
the cold stellar material stalls, as a consequence of the drop- 
ping cocoon pressure. This is equivalent to assuming a con- 
stant volume of the cocoon cavity inside the star. The pres- 
sure drop for the relativistic cocoon can be derived through 
dE c = —e c S e c s dt where E c is the area of the free surface 
through which the cocoon material expands, e c the cocoon 
energy density and c s = c/ a/3 is the sound speed of the rela- 
tivistic gas. Writing the cocoon volume as V c ~ £ e r*, we can 
obtain the pressure evolution: 



: Pc.br exp 



ct 



V3r, 



(8) 



where p c b r is the cocoon pressure at the moment of shock 
breakout. 

As the cocoon pressure decreases, fresh jet material pass- 
ing through the cocoon is less tightly collimated. Under isen- 
tropic conditions the jet Lorentz factor increases linearly with 



= Oj, br exp 



ct 



4%/3r, 



(9) 



Dissipative jet propagation gives similar results (with merely 
a different numerical coefficient ^ 0(1) inside the exponen- 
tial), provided that V varies roughly as a power of £_,-. The 
angular distribution of integrated energy, as observed in the 
afterglow phase, is computed by integrating the instantaneous 
luminosity per unit solid angle from the moment the jet be- 
comes visible along a given line of sight (?i. . s .) until the end 
of the burst: 



dE 
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dh 
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(10) 



fi.o.s. is obtained by inverting eq. [9] Such integration is 
valid provided that the jet opening angle at time 7grb is 
smaller than the natural opening angle of the jet: 7grb < 
4v / 3('v/c)log(#o/6'br)- For the fiducial numbers assumed 
(r* = 10" cm, 0y,br =1° an d #o = 30°) this corresponds to 
~ 100 s comoving burst duration. Interesting effects should 
be expected for longer bursts, as discussed in § 4. Assuming 
a jet with constant luminosity L and for all the lines of sight 
that satisfy fbr < fi.o.s. <^ ?grb, eq.[l0]gives the jet structure 
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and dE/dil ^constant inside the core radius 6,-^. This an- 
gular dependence, which characterizes a "structured jet" or 
"universal jet" (Rossi et al. 2002, 2004; Zhang & Meszaros 
2002; Salmonson 2003; Lamb, Donaghy & Graziani 2003), is 
of high theoretical interest. Jets with this beam pattern repro- 
duce afterglow observations. If the jet is powered by fall-back 
of material from the star to the accretion disk, the mass accre- 
tion rate would be anti-correlated with the radius of the star, 
for a given stellar mass. This would produce a roughly con- 
stant value of Ljr+ that would reproduce the so-called Frail 
relation (Frail et al. 2001; Panaitescu & Kumar 2001) as a 
purely viewing-angle phenomenon. In this case all jets would 
be alike, but different observers would see them from different 
angles, deriving different energetics. 

In the above equations we have assumed for simplicity that 
the jet reaching the surface of the star is uniform. As shown 
by simulations (e.g. Zhang et al. 2003), it is more likely that a 
Gaussian jet emerges from the star. On the contrary, boundary 
layers may be produced by the interaction of the jet with the 
collimating star, resulting in edge brightened jets (see the in- 
set of Fig.^. It can be shown easily that the 8~ 2 pattern does 
not depend on the assumption of uniformity. Fig.^shows the 
integrated energy distribution for uniform, Gaussian and hol- 
low intrinsic jets. Even though small differences are present at 
the edges (the jet core and the outskirts), the general behavior 
is always dE/dfl oc 6~ 2 . 

4. DISCUSSION 

Our prediction that the jet opening angle evolves with time 
has two major, in principle testable, consequences. First, the 

2 Note that, if the jet is causally connected at breakout as suggested by 
eq.|6| the jet would freely expand to an angle 6j = 1 /Vj\x > dj.bi of the order 
of a few degrees. This would result in an initially co nsta nt opening angle. 
The only effect on the final energy distribution of eq.ITTlis to increase the 
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brightness of a typical GRB should tend to decrease with 
time, assuming a constant efficiency of gamma-ray produc- 
tion, since the photon flux at earth is proportional to the en- 
ergy per unit solid angle. This behavior is certainly seen in 
FRED (Fast Rise Exponential Decay) single pulsed events, 
while it is harder to make a quantitative comparison in the 
case of complex light curves. Some events, like the bright 
GRB 990123, have multi-peaked lightcurves which show a 
clear fading at late times. On the other hand there exist rare 
events, such as GRB 980923, in which an increase of lumi- 
nosity with time is observed. Lacking afterglow observations, 
it is hard to determine whether such events are peculiar. 

To understand the second consequence of the model, con- 
sider an observer at an angle O \, S > &bi where 6*b r is the jet 
opening angle at breakout. Initially, this observer lies out- 
side the beaming cone and does not detect the GRB. Later, 
as the jet spreads beyond 6> b s , the observer detects the burst. 
The beginning of the GRB emission is therefore observer de- 
pendent. This causes a correlation between the burst duration 
and energetics: the longer the burst the larger its detected en- 
ergy output. A correlation between burst duration and fluence 
is detected in the BATSE sample but whether it is intrinsic 
or due to systematic effects is still a matter of debate. Fi- 
nally, if any isotropic emission should mark the jet breakout 
(MacFadyen & Woosley 1999; Ramirez-Ruiz, MacFadyen & 
Lazzati 2002), such emission would not be followed imme- 
diately by 7-rays for observers lying outside of 0^. This de- 
lay may explain the unusually long delays between precursors 
and main emission detected in several BATSE GRBs (Lazzati 
2005). 

Some complications can arise for very compact stars and/or 
very long-lived jets. In computing the energy profile (~eq. ll It 
we have implicitly assumed that the jet expands exponentially 
until it dies. This approximation obviously fails if the engine 
is long-lived, since at some point the jet reaches its natural 
opening angle - determined, e.g., by the accretion disk geom- 
etry - and the expansion stops. Alternatively for some com- 
bination of the parameters, the cocoon may occupy a small 
opening angle (eq. |5j . As the jet expands, it eventually hits 
the cold star material and a different expansion law, driven 
by the jet pressure onto the stellar material, sets in. Again, 
the expansion of the jet is slowed down and as a consequence 
a shallower distribution of energy with off-axis angle is ex- 
pected. Such modifications to the 0~ 2 law would manifest 
themselves in the late stages of the afterglow evolution. The 
higher than expected energy at large angles would flatten the 
light curve decay or produce a bump in the afterglow. Such 
behavior has been claimed in several radio afterglows (Frail et 
al. 2004; Panaitescu & Kumar 2004), and could also explain 

3 Note that a simple two-component jet cannot explain the fast variability 
observed in the optical lightcurve. It can only reproduce the overall behavior 



the two breaks observed in the light curve of GRB 030329 3 . 
If the outflow reaches its natural opening angle or crashes into 
the cold stellar material while still active, the time-integrated 
jet can be described as the sum of a 9~ 2 jet plus a roughly 
uniform jet with a larger opening angle — in other words, a 
two-component jet (Berger et al. 2003). This behavior is tan- 
talizing also if we consider the interaction from the point of 
view of the star. The expansion of the cocoon drives a shock 
wave into the cold stellar material where it deposits approxi- 
mately 10 51 erg (Zhang et al. 2003; Lazzati & Begelman in 
preparation). If the cocoon is small and the jet comes into 
contact with the star while still active and expanding, addi- 
tional energy is given to the star. Complex afterglows should 
therefore be associated with particularly energetic hypernovse, 
such as SN2003dh. 

5. SUMMARY AND CONCLUSIONS 

In this Letter we have analyzed the evolution of a GRB jet 
during its transition from stellar confinement to free expan- 
sion. We focused mainly on the resulting angular structure 
of the outflow, and reached two main conclusions. First, the 
angular structure of the jet in the afterglow phase does not 
necessarily reflect the way in which the energy is released by 
the inner engine. If the jet opening angle is not constant dur- 
ing the ~ 100 s of the engine activity, the angular structure of 
the jet results from the integration of this evolution and has 
little to do with the pristine beam pattern. Second, we ana- 
lyzed the most probable evolution in the collapsar scenario. 
We find that the interaction of the jet with the surrounding co- 
coon causes the opening angle to increase exponentially with 
time. This evolution, for a constant energy release from the 
inner engine and independently of the intrinsic beam pattern, 
produces an angular structure with dE/dfl oc 9~ 2 . Such an 
energy distribution produces a broken power-law afterglow 
(Rossi et al. 2002, 2004; Salmonson 2003, Kumar & Gra- 
not 2003; Granot & Kumar 2003). It can also explain the so- 
called Frail relation (Frail et al. 2001; Panaitescu & Kumar 
2001) if Ljr^/c is roughly constant for different GRBs. The 
origin of this energy distribution is here explained for the first 
time in the context of hydrodynamical fireballs (see Lyutikov 
et al. 2003 for a discussion on the context of electromagnet- 
ically dominated fireballs). We briefly discuss in § 4 some 
testable prediction of this model. 

We thank Miguel Aloy, Gabriele Ghisellini and Andrew 
MacFadyen for useful discussions. This work was supported 
in part by NSF grant AST-0307502 and NASA Astrophysical 
Theory Grant NAG5- 1 2035 . 

of the optical and radio lightcurves. 
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